Simulation of vector fields with arbitrary second-order correlations.

Temporally-stationary electromagnetic fields with arbitrary second-order coherence functions are simulated using standard statistical tools. In cases where the coherence function takes a commonly-used separable form, a computationally-efficient variation of the approach can be applied. This work provides a generalization of previous spatio-temporal simulators which model only scalar fields and require either restrictions on the coherence function or consider only two points in space. The simulation of a partially-polarized Gaussian Schell-model beam and a partially-radially-polarized beam are demonstrated.